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Abstract. Proton-proton and proton-ry invariant mass distributions and the total cross section for the pp — > 
pprj reaction have been determined near the threshold at an excess energy of Q = 10 MeV. The experiment 
has been conducted using the COSY-11 detector setup and the cooler synchrotron COSY. The determined 
invariant mass spectra reveal significant enhancements in the region of low proton-proton relative momenta, 
similarly as observed previously at higher excess energies of Q = 15.5 MeV and Q = 40 MeV. 
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1 Introduction 

The complexity of the structure of hadrons constitutes 
the basic difficulty in the quantitative description of the 
hadronic interaction in the medium energy regime. There- 
fore, this interaction is not well understood especially in 
the meson-nucleon and meson-meson sector, where an ad- 
ditional difficulty is the relatively poor experimental data- 
base. Particularly challenging are investigations of inter- 
actions involving flavour-neutral mesons. This is due to 
the short life-time of these mesons which can neither be 
used as targets nor as beams. Thus, in practice such in- 
teractions can be accessed only indirectly via observables 
like excitation functions or invariant mass distributions. 
Measurements of these observables are especially useful in 
the close-to-threshold region where the final state parti- 
cles are produced with low relative velocities. Among the 
basic flavour neutral mesons the rj is of particular interest 
since its interaction with nucleons is strong enough to be 
detectable with the presently achievable experimental pre- 
cision p], and since its interaction seems to be sufficently 
strong to form an eta-mesic nucleus 01!]. The existence 
of such kind of nuclear matter is vividly discussed [4] and 
there are ongoing experimental programs searching for a 
signal of such a state [5]. 

The earlier high statistics measurements of the pp — > 
pprj reaction at an excess energy of Q = 15.5 MeV from 
the COSY-11 collaboration [5], and also the measurements 
of the TOF group [7] at Q = 15 and 41 MeV, revealed 



that there exist significant enhancements in the invari- 
ant mass distributions of pp and prj subsystems at higher 
values of proton-proton invariant mass and lower values of 
the proton- rj invariant mass. One of the plausible explana- 
tions for these enhancements could be an influence of the 
proton- rj interaction [T][8] . If this is the case one could use 
such observables for the estimation of the strength of this 
interaction. However, the observed invariant mass distri- 
butions could be also plausibly explained by contributions 
of higher partial waves [UlllOj or by an energy dependence 
of the primary production amplitude [llill2j . Therefore, 
in order to verify the correctness of the proposed explana- 
tions it is of importance to investigate the dependence of 
the strength of the enhancements as a function of the ex- 
cess energy. Qualitatively, with decreasing excess energy 
the contribution from the higher partial waves should de- 
crease whereas the influence of the interaction should be 
more pronounced. 

Certainly, most effectively, contributions from higher 
partial waves could be disentangled by the determination 
of the analysing powers and spin transition coefficients |10[ 
I13U14] , yet such investigations are not planned in the near 
future at COSY which is at present the only laboratory 
where it can be conducted. This makes the determination 
of the energy dependence of the distributions of the differ- 
ential cross section for the pp — > pprj reaction even more 
important for studies of the proton-77 hadronic interaction 
and for studies of the properties of nucleon resonances [HI 

eng. 
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In this article we present distributions of the proton- 
proton and proton-ry invariant masses at the excess energy 
of Q — 10 MeV which is significantly closer to the thresh- 
old with respect to the previous studies. Although the 
original experiment at Q = 10 MeV has been devoted to 
the investigations of the analysing power for the pp — > ppr\ 
reaction [16 and has been performed with a polarised pro- 
ton beam, the data enable also the determination of the 
spin averaged observables after appropriate offline " depo- 
larisation" of the beam explained in section 12.11 In Sec- 
tion[2]we briefly describe the experimental set-up, present 
the experimental principles of the measurement, and de- 
scribe the method of the data analysis. In section [3] the 
determined spectra are compared to the analogoues re- 
sults determined at the excess energy of Q = 15.5 MeV, 
and the final conclusions are drawn. 



2 Experiment 

The measurement of the pp — > pprj reaction has been per- 
formed at the cooler synchrotron COSY [17] at the Re- 
search Center Julich in Germany using the COSY-11 de- 
tector setup [18] , presented schematically in Figure [1] 




Fig. 1. Schematic view of the COSY-11 detector setup [18] . 
Dl and D2 represent drift chambers. SI, S3, S4 denote scintil- 
lator counters. Simon is the silicon strip detector used for the 
detection of the elastically scattered protons. Superimposed 
solid lines indicate final state protons from the pp — > ppr\ reac- 
tion. The size of detectors and their relative distances are not 
to scale. 



The proton beam with a momentum of 2.010 GeV/c, 
corresponding to an excess energy of Q = 10 MeV, has 
been scattered on H2 molecules from an internal cluster 
jet target [I9 ] l20 j . installed in front of the COSY magnet. 
Reaction products carry lower momenta than the beam 
protons, therefore are bent more in the magnetic field of 



the dipole magnet. Positively charged ejectiles leave the 
scattering chamber through a thin exit window reaching 
the detection system operating under atmospheric pres- 
sure. The hardware trigger was based on signals from 
scintillator detectors only. It was adjusted to register all 
events with at least two positively charged particles. For 
this aim coincident signals in the SI and S3 detectors were 
required. In the case of the SI detector only these events 
were accepted for which either two separate segments were 
hitted or an amplitude of the signal in a single module was 
higher than the certain threshold value. Based on the data 
analysis from the previous experiments, the threshold was 
set high enough to reduce significantly the number of sin- 
gle particle events, and at the same time to accept most 
events (almost 100%) with two protons passing through 
one segment. Next, in the off-line analysis it was required 
that at least two tracks are reconstructed from signals 
measured by means of two planar drift chambers Dl and 
D2. The trajectories of the positively charged particles re- 
constructed in drift chambers are further traced through 
the magnetic field of the dipole back to the interaction 
point. In this way the momenta of the particles can be 
reconstructed with a precision of 6 MeV/c (standard de- 
viation) [6]. The timc-of- flight measurement between the 
scintillator hodoscope SI and the scintillator wall S3, and 
the independently reconstructed momentum enable a par- 
ticle identification by means of the invariant mass tech- 
nique. The COSY-11 mass resolution allows for a clear 
seperation of groups of events with two protons, two pions, 
proton and pion and also deutron and pion pQ . Further on 
the produced meson is identified using the missing mass 
method. A more detailed description of the method and 
results of the identification of the pp — > pprj reaction can 
be found in reference 



2.1 Off-line depolarisation of the beam 

Originally the experiment was dedicated to the measure- 
ment of the beam analysing power for the pp — > pprj re- 
action [16, 21 . Therefore, the proton beam has been ver- 
tically polarised. The vertical polarisation of the beam is 
defined as an asymmetry of populations of particles in the 
spin up (N+) and down (AT_) states with respect to the 
vertical axis, integrated over the whole period of measure- 
ment: 

N + - Y_ 



Py = 



N+ + N. 



(1) 



In the discussed experiment the direction of the polarisa- 
tion was being flipped from cycle to cycle. Hence, for the 
so called "spin up cycles" we define the spin up polarisa- 
tion as: 



FT 



L +,i 



Ei=«p + Ei 



(2) 



and analogously for " spin down cycles" the spin down po- 
larisation reads: 



pi 



dn ' 
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(3) 
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where, n+j and ri— i denote the number of protons in the 



:th 



cycle, in spin up and down state, respectively. Please 
note, that according to above definitions the following re- 
lations are valid: 



i—up 
i—up 



i—dn 

£ 

i—dn 



N— 



(4) 



the geometry of the COSY-11 detector setup including the 
conditions of the beam and target and the resolution 
and efficiency of the detectors pQ. The simulated events 
were analysed with the same program as the experimen- 
tal data. Subsequently, functions of the type: 



/(mm 2 ) = a ■ fpp-,pp2^{mm 2 ) + f3 ■ fp P ^ P p3iv(mm 2 ) 
+ 7 ' ipp—ppi-nimm 2 ) + S ■ fpp-,p m {mm 2 )i 



(7) 



The beam can be effectively depolarized e.g. by assigning 
to the events in spin up cycles the weights w which can 
be derived from the requirement that the numerator of 
Equation [T] has to vanish: 



i—dn 



i—dn 



Thus, combining Equations [2] and [3] with Equation [5] we 
obtain the following formula for the value of w: 



p 1 E, 



dn ■ 



up ' 



(6) 

where L re i := jj denotes the relative luminosity for the 
spin up and down cycles. Taking into account the numer- 
ical values of pT = 0.658 ±0.008, P l = 0.702 ± 0.008, and 
L rei = 0.98468 ± 0.00056 [l6pT] one gets w = 1.083. 

The weight w, assigned to events in spin up cycles, 
does not change the absolute value of cross sections, as 
the same weight has been applied in both: the calculation 
of the number of events originating from the pp — ► pprj 
reaction and the determination of the luminosity from the 
pp — > pp elastic scattering. 



2.2 Data analysis 

The method applied to the determination of the differen- 
tial cross sections follows the procedures described in [6], 
therefore for any details the reader is reffered to that pa- 
per. Here we shall only briefly describe the main steps of 
the data analysis and emphasize the differences between 
the methods used in both studies. 

After particle identification we continued the analy- 
sis with the depolarization of the experimental data, ac- 
cording to the procedure described in Section 12.11 Next, 
we determined the covariance matrix and performed the 
kinematical fitting [TJG)]. 

For the description of the relative motion of the pro- 
tons and the rj meson, following reference [6], we have 
chosen the squares of the invariant masses - s pp and s pri - 
of the proton-proton and proton-77 systems, respectively. 
Optimizing the statistics we have divided the range of s pp 
and s p „ into 20 bins. For each bin of these variables we 
have determined the spectrum of the square of the missing 
mass. Analogoues spectra were simulated for the pp — > ppr\ 
reaction and for the background channels. The simulation 
program, based on the GEANT3 [22] code, accounts for 



were fitted to the data, with a, (3, 7, and S treated as free 
parameters responsible for the normalisation of the sim- 
ulated missing mass spectra (/) of reactions indicated in 
subscripts. In order to determine the background free in- 
variant mass distributions the experimental missing mass 
squared spectra were fitted separately for each bin of s pp 
and Sp V . As an example, the missing mass distributions 
for arbitrarily chosen bins of s pp and s p „ are presented in 
Figure [2] 
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Fig. 2. Missing mass squared distributions for arbitrarily cho- 
sen bins of Spp (left) and s pv (right), as measured at the excess 
energy of Q = 10 MeV. Dots represent the experimental data 
points along with their statistical erros, whereas the solid line 
is the best fit of the sum of the signal and background as ob- 
tained in the Monte-Carlo simulations. The shaded part shows 
the generated multi-pionic background. 
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The mumbers of r\ mesons N v for the individual inter- 
have been calculated as: 



vals of s pp and s 



at the excess energy of Q 
right part of Figure [3] 



15.5 MeV [6], displayed in the 



N„ = 5 ■ 



fpp^pp-nimm 2 ) d(mm 2 ), (8) 



and the statistical errors cr(N v ) have been estimated as: 



a(N v ) = a(6) 



fpp^pp V (mm 2 ) d(mm 2 



(9) 



where cr(<5) are the estimates of the S parameter uncer- 
tainties (standard deviations) determined by means of the 
MINUIT minimization package [24] . 

The systematic error of N v has been estimated to be 
not larger than 8%, based on the dependence of the re- 
sults on different assumptions for i) background estima- 
tion (~2%), ii) description of the proton-proton final state 
interaction (~5%), and iii) inaccuracy in efficiency for re- 
construction of both proton trajectories (~6%). The un- 
certainty in the number of the background events under 
the peak of the 77 meson was estimated as differences be- 
tween results obtained by fitting to the background a) the 
first order polynomials, b) the second order polynomials, 
c) the sum of two Gaussian functions, and d) the distri- 
butions simulated for the multi-pion production [31~1[2"T] . 
The inaccuracy due to the model used for the description 
of the proton-proton FSI was estimated conservatively as 
a difference in results determined when using parameter- 
ization of the proton-proton S-wave interaction [25j and 
when neglecting the FSI and taking into account a homo- 
geneous phase-space distribution of the momenta of final 
state particles. 

The determination of the luminosity was based on the 
comparison of the measured differential pp — > pp elas- 
tic scattering cross sections to the data of the EDDA 
group [26] . For the detailed method of the luminosity and 
acceptance calculation the interested reader is referred 
to [T]l23j. The integrated luminosity was extracted to be 
L = 58.53 ±0.03 nb" 1 . 



3 Results and conclusions 

The total cross section evaluated as an integral of the s pp 
distribution equals to a — 1.27±0.04±0.13 /jb, where the 
first error is the statistical and the second the systematic 
one. The latter accounts for the quadratic sum of inde- 
pendent contributions from 8% systematic error of N v , 
3% systematic error of the luminosity determination [6], 
and 6 % uncertainty in the acceptance estimation 6J . The 
determined total cross section is in line with the previous 
measurements performed independently by various exper- 
imental groups [6"ll2"T]. 

The results on the differential cross sections for the 
PP ~ * PPV reaction as a function of the square of the 
proton-proton and proton-77 invariant masses are given in 
Table [1] and presented in the left part of Figure [3] They 
are compared with the differential cross sections measured 






- phase space 
FSI„ 




■ -H 



3.53 3.54 3.55 3.56 

s BD [GeVV] 



phase space 
FSI nr 



3.54 3.5f 
s dd [GeV 2 /c 4 ] 




phase space 




fsi --X 


\ 






/ 0' 


\ * 

\ k 


Av 




{*/ 





2.22 2.24 

s m [GeV 2 /c 4 ] 



Fig. 3. Distributions of the square of the proton-proton and 
proton-77 invariant masses as measured at Q = 10 MeV (left) 
and Q = 15.5 MeV (right). Solid lines represent the homo- 
geneous phase-space distributions, while the dotted lines are 
the theoretical predictions taking into account the 1 So proton- 
proton final state interaction. 



The homogeneous phase-space distributions, shown in 
Figure [3] by solid lines, completely disagree with the ex- 
perimental data for all distributions presented. The peak 
observed at small values of s pp is associated with strong 
proton-proton final state interaction (FSI). On the other 
hand the dotted lines, which represent the phase-space dis- 
tribution convoluted with the proton-proton FSI describe 
the data quite well in the range of small invariant masses 
of the proton-proton system and in the range of large s pv , 
but for both excess energies there is a significant deviation 
of the theoretical predictions from the experimental data 
in the range of large values of s pp and small values of s pv . 
In the calculations we have used the parameterization of 
the proton-proton FSI given in reference [28] . The curves 
including the proton-proton FSI have been arbitrarily nor- 
malized in the range of small values of s pp and close to the 
upper limit of the s pv distribution. 

A preliminary result from a comparative analysis of 
the pprj and ppr( system indicates that the observed en- 
hancement is rather not due to the meson-proton interac- 
tion [3Tfll3"T] . Preliminary results show that the enhance- 
ments in the invariant mass distributions are also present 
in case of the pp — > pprj' reaction [30]l3T] . Due to the fact 
that the interaction between the rj' meson and proton is 
more than an order of magnitude weaker than the one be- 
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Table 1. Distribution of the square of the invariant mass 
of the proton-proton and proton-77 systems for the pp — > ppr\ 
reaction at the excess energy of Q = 10 MeV. 
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Another explanation for the bumps observed at higher 
values of s pp were put forward by DelofT [11 who explains 
the observed spectra by allowing a linear energy depen- 
dence of the leading 3 Pq -^ 1 SqS partial wave amplitude. 
Recently also Ceci, Svarc and Zauner [T2"l[3"2"] have shown 
that the excitation function and the enhancement in the 
invariant mass spectra can be very well described by the 
energy dependence of the production amplitude when the 
negative interference between the 7T and the 77 meson ex- 
change amplitudes is assumed. However, for the quanti- 
tative confirmation of these hypotheses still more precise 
data on the energy dependence of the enhancement in the 
invariant mass spectra are required. 
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